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Abstract 

While first-order channel statistics, such as bit-error rate (BER) and outage probability, play an important role in the design 
of wireless communication systems, they provide information only on the static behavior of fading channels. On the other hand, 
second-order statistics, such as level crossing rate (LCR) and average outage duration (AOD), capture the correlation properties of 
fading channels, hence, are used in system design notably in packet-based transmission systems. In this paper, exact closed-form 
expressions are derived for the LCR and AOD of the signal at a receiver where maximal-ratio combining (MRC) is deployed over 
flat Rayleigh fading channels in the presence of additive white Gaussian noise (AWGN) and co-channel interferers with unequal 
transmitting powers and unequal speeds. Moreover, in order to gain insight on the LCR behavior, a simplified approximate 
expression for the LCR is presented. As an application of LCR in system designs, the packet error rate (PER) is evaluated through 
finite state Markov chain (FSMC) model. Finally, as another application again by using the FSMC model, the optimum packet 
length to maximize the throughput of the system with stop-and-wait automatic repeat request (SW-ARQ) protocol is derived. 
Simulation results validating the presented expressions are provided. 


I. Introduction 

The first order statistics are commonly used as important performance metrics, however they lack information for the overall 
system design and performance since they only capture the static behavior of the channel. Outage probability, the probability 
that the received signal-to-interference plus noise ratio (SINR) is below a certain threshold, and bit error rate (BER), the rate 
at which errors occur during the transmission, are such static metrics. On the other hand, the second-order statistics reflect the 
correlation properties of the fading channel providing a dynamic representation of the system performance. The level crossing 
rate (LCR) and the average outage duration (AOD) are examples of second-order statistical metrics and have been used in 
many applications for designing and evaluating the performances of wireless communication systems. They have been used in 
finite-state Markov chain (FSMC) channel modeling ||T), the analysis of handoff algorithms Q, Markov decision process (MDP) 
and partially observable Markov decision process (POMDP) formulations 0-1^, and packet error rate (PER) evaluation |[^. 
Moreover, since AOD determines the average length of error bursts, it plays an important role in choosing the packet length 
for uncoded and coded systems, designing interleaved and non-interleaved coding schemes ||7), determining the buffer size and 
transmission rate for adaptive modulation schemes Q, and estimating the throughput of different communication protocols 
such as automatic repeat request (ARQ) |[^. Since LCR and AOD are closely related metrics, it is imperative to obtain the 
exact closed-form expressions for both LCR and AOD. 

In his pioneering work. Rice obtained the LCR expression for a Rician-distributed fading signal in a single-user additive 
white Gaussian noise (AWGN) channel ng- Many works have been conducted later on to find LCR and AOD expressions in 
noise-limited systems covering different combining schemes with correlated and uncorrelated branches HD-®- Fewer works 
addressed interference-limited systems due to its relative complexity and the hardship of finding the joint probability density 
function (PDF) of the signal-to-interference ratio (SIR) and its time-derivative, and its challenging integration steps. Examples 
of such works include E) which derived the LCR and AOD expressions for interference-limited selection combining (SC) 
system over Rayleigh distributed channel while the interferers’ channels are Rician distributed. In | [T8) , Yang and Alouini 
derived the LCR expression where the desired signal is received over i.i.d. diversity paths and are Nakagami, Rician or 
Rayleigh distributed and so are the interferers; however, background noise effect was ignored. Recently, Fukawa et al. Q 
derived an approximate LCR expression for a point-to-point system subject to frequency-selective Rayleigh fading with the 
receiver deploying maximal-ratio combining (MRC) and the interfering signal components are of unequal powers. Seldom 
attention was given to systems where noise and interference powers are comparable. In E) an exact closed-form expression 
was obtained for the LCR of a single-antenna receiver with multiple interferers traveling at different speeds with unequal 
powers under AWGN and Rayleigh fading channel. However, it cannot be applied to MRC and hence loses the diversity gain. 
In this paper, the exact closed-form expressions of LCR and AOD are derived for a more general setting where the receiver 
with multiple antennas uses MRC diversity combining in the presence of multiple co-channel interferers that have different 


powers and speeds while the effects of both AWGN and the Rayleigh flat fading channels are taken into consideration. To the 
best of our knowledge, this is the first time to report these exact closed-form LCR and AOD expressions for this setting. It is 
worth noting that the results in |[^ are the special cases of the expressions derived in this paper when the number of receive 
diversity branches is set to 1 as shown in Section |V-C A summary of the most relevant works in the literature on the LCR 
derivation is listed in Table U 


TABLE I: A summary of the scenarios solved w.r.t. LCR derivation. 


Problem 

Contribution 

Research Papers 

Single user Rician fading channel 

Exact closed-form solution (ECFS) 

1101 

Correlated Rayleigh-fading multiple-branch equal-gain 
predetection diversity combiner 

Approximate closed-form solution (ACES) 

HI] 

Independent Rayleigh-fading with 



two-branch selection combining (SC), equal-gain (EGC) combining 
and maximal-ratio combining (MRC) for a single-user AWGN channel 

ECFS 

m 

Independent Nakagami-fading multiple-branch SC, EGC and MRC 

ECFS for SC and MRC 

Numerical integration solution for EGC 

(H) 

EGC and MRC over independent multiple-branch 

ECFS for MRC with Rayleigh i.i.d. branches 

IB 

with different fading distributions 

Numerical integration solution for the other scenaiios 

Independent unequal-power multiple-branch 
Rayleigh-fading MRC 

ECFS 

IB 

Independent unequal-power multiple-branch 

Weibull-fading SC 

ECFS 

IB 

Interference-limited SC system with Rayleigh fading desired user 


Ozl 

and i.i.d. Rician-distributed co-channel interferers 

ECFS 

with equal average powers 


Interference-limited MRC system with i.i.d. Nakagami, 

Rician or Rayleigh distributed desired signal branches 

ECFS 


and i.i.d. co-channel interferers 



Interference-limited MRC point-to-point system with i.i.d. 
frequency-selective Rayleigh fading channel with 

ECFS for equal average power case 

0 

equal and unequal average power interferers 

ACFS for unequal average power case 

Single antenna nodes with Rayleigh-fading channels with interferers 
having unequal average powers and different speeds 

ECFS for the LCR of the SINR 



LCR and AOD take place in a variety of applications ranging from estimating the system PER by using the ESMC model 
of the system to optimizing the packet length of automatic repeat request (ARQ)-based systems. In pO) , Yousehzadeh and 
Jafarkhani estimated the PER of a multi-input multi-output (MIMO) system in the presence of co-channel interferers under 
AWGN and fading channels by using ESMC model although strong approximations were made in the SINR expression, while 
Eukawa et al. @ used the ESMC model to evaluate the PER of the point-to-point system under frequency selective fading. 
Many works have been done in literature on automatic repeat request (ARQ)-based systems due to their practical importance 
in improving the data transmission efficiency ED- Different ARQ schemes, including stop and wait (SW) ARQ, go-back-N 
(GBN) ARQ, selective-repeat (SR) ARQ, and hybrid (H) ARQ, have been analyzed for different point-to-point AWGN and 
fading channels Q, ||22)-||2^. In Q, the authors evaluate the throughput using ESMC modeling of the system under GBN 
and SR ARQ protocols. The system’s throughput is analyzed for SW ARQ-based system under slow Rayleigh fading with 
adaptive packet length in p2) , while it is analyzed for GBN ARQ-based system under Rician fading using Gilbert-Elliott 
channel (GEC) model in |23| . The throughput is analyzed and maximized w.r.t. information rate for HARQ-based system 
over an orthogonal space-time block coding (OSTBC) MIMO block fading channel in p5) . ESMC is used in | [24) to model 
the packet error process in an SR ARQ-based system under Rayleigh fading, and the delay statistics is investigated. While 
those works considered single-user systems for which ARQ was originally developed, recently more attention was given to 
multi-user (ML!) systems deploying ARQ protocols such as fZT) which considers linear precoder optimization to maximize the 
throughput in a MU-MIMO HARQ-based system while p8) deals with the user scheduling aspect to improve the throughput. 
p9| considers minimizing the transmission power via resource allocation in LTE uplink where HARQ is utilized assuming 
block fading Rayleigh channel model. It is worth mentioning that these recent works, p7)-p9), do not consider the channel 
time-correlation aspect. 

Since the throughput is an important metric in communications systems, many works tackled the throughput maximization 
problem in different settings under different parameters, e.g. [30) maximizes the throughput of an abstract system deploying 
one-bit ARQ suffering from delay through multi-user scheduling where POMDP framework is adopted. An important parameter 
affecting the throughput, specially in ARQ systems, is the packet length. Packet length and rate adaptation is tackled in pT[ to 
maximize the throughput in wireless LANs assuming Nakagami-m fading while no retransmissions are allowed, hence, allowing 
constrained packets loss. |3^ proposed an algorithm to maximize the throughput in Bluetooth piconets through selecting the 






















optimum packet length given a finite selection set. The authors in p3| maximize the throughput of a point-to-point system 
utilizing truncated ARQ scheme under flat Rayleigh fading, where the optimum packet length and modulation size are chosen 
out of a given finite set based on an iterative suboptimum algorithm. In again a point-to-point uncoded transmission 
system utilizing SW-ARQ with infinite retransmissions is studied under slow Rayleigh flat fading and an expression for the 
optimum data packet length maximizing the throughput was derived though the results are restricted to systems where the 
BER has exponential form. Again, we note that these works do not consider the channel time-correlation in the analysis. 

In this paper, the exact LCR and AOD expressions are derived, in contrast to e.g., 1^, under the effects of both AWGN 
and fading channels, in contrast SIR or signal-to-noise ratio (SNR) focused works e.g., | |14| , HI). with a receiver employing 
MRC, in contrast to GH- Further, to easily acquire the LCR behavior with lesser computational complexity under different 
conditions, an approximate LCR expression is derived. By developing an FSMC model that makes use of the exact closed-form 
LCR expression and captures the channel time-correlation via incorporating the LCR in the model, the PER of the system, 
and the optimum data packet length for throughput maximization for a multi-user SW-ARQ system are derived. The main 
contributions of this work are summarized as follows. 

• Exact closed-form expressions for the LCR and AOD of the received SINR for an MRC receiver are derived in the 
presence of AWGN and multiple co-channel interferers with unequal powers, and moving with unequal speeds subject to 
Rayleigh flat fading channels. 

• Approximate closed-form expressions of LCR and AOD are derived. These simplified expressions can help system 
designers streamline the design process without explicitly evaluating the LCR with every change of system parameters. 
Moreover, the computation complexity of the LCR is significantly reduced, especially for high number of interferers and 
receive antennas. 

• The system is modelled by FSMC via the exact LCR expression and closed-form PER is derived. 

• The formula for optimum packet length maximizing the throughput of the SW-ARQ system is derived. 

Finding an exact closed-form solution for the LCR is a difficult process in general due to involved derivations including the 
joint PDF of the SINR and its time derivative along with difficult integrals that are challenging to be transformed to easier 
integrals. The difficulties arise from the following aspects included in this paper. 

• Both noise and interference are considered. As well known, noise-limited systems are rather easier since only SNR 
is tackled and hence the interferers’ channels in the SINR denominator are absent. Again as well known, although 
interference-limited systems are more difficult than noise-limited systems, it is still an easier problem than the one 
incorporating both noise and interference. This is evident from the abundance of works considering either noise or 
interference and scarcity of works considering both. 

• The interferers are assumed to be having different transmitting powers and speeds that lead to mathematical difficulties 
in obtaining the joint PDFs and in solving the integrations. 

• An integral emanate corresponding to each interferer. Since a solution for any number of interferers is expected, an 
analytical solution for arbitrary number of integrations is not trivial contrary to the approaches that use numerical solutions 
for integrations. 

• The receiver has multiple antennas and deploys MRC which constitute complicated integrations. 

The theoretical expressions are validated by comparisons with the simulation results under varying system parameters. Some 
interesting findings in this work are 

• Given a fixed ratio of the SINR threshold to the average SINR, i.e. normalized threshold 7th/7avg^ the LCR value of the 


received SINR is independent of the value of the desired signal’s power relative to the interferers’. Section VII-A 


The exact LCR expression is very accurate in terms of matching the simulation results. Section VII 


• The PER obtained from the LCR expression via the FSMC modeling is very accurate in terms of matching the simulation 
results. Section [Ix| 

The paper is organized as follows. The system model is presented in Section The joint PDF of the SINR and its time 
derivative are derived in Section Iffl and used to obtain the LCR in Section IV The LCR expressions of simplified systems are 
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derived as special cases in Section V^ and are compared with previously reported expressions in the literature. The simplified 
approximate LCR expression is derived in Section VI Section |VII| presents the validation of the derived expressions via 
numerical results along with the investigation of the system parameters’ effects on the LCR. The AOD derivation is given in 
Section VIII and verified via numerical examples. As an application of the new exact LCR expression, the system is modeled 
as an FSMC model and the PER is evaluated in Section IX Another application of the LCR expression is presented in Section 
jXj where the throughput of the system under SW-ARQ protocol is maximized by obtaining the optimal packet length as a 
function of the system parameters which includes the LCR. Finally, conclusions are drawn in Section [^ 


II. System Model 

We consider a mobile radio system with N + 1 transmitter nodes and a receiver node. Each transmitter and the receiver 
have a single antenna and L antennas, respectively. N transmitter nodes are interfering users, thus there is a single desired 
transmitter. All channel paths in the system are subject to flat Rayleigh fading with unit average power. The time autocorrelation 









of each channel path is given by Jo{27rfnT), where Jo{.), fo, and t are the Bessel function of the zeroth order, the maximum 
Doppler frequency associated with this path and the time difference between the two correlated samples, respectively | [^ . 
The desired transmitter power and the interfering transmitter power are and p„, Vn G Af = {1,..., iV}, respectively. The 
receiver performs MRC over the received signal. The received noise z is an additive white Gaussian noise (AWGN) with zero 
mean and covariance matrix i?{zz^} = Nol. This model can be generalized to any system with L i.i.d. diversity paths, where 
again the receiver performs MRC over the L diversity paths. The diversity can be also achieved in time or frequency, not 
necessarily in space. 

Let hD(f) and h„(f) denote the L x 1 channel vectors between the receiver and the desired user and ^th 

interfere!', respectively. 

The received signal before applying MRC, y{t), is given by 


N 

y{t) = hD{t)so{t) + ^ K{t)sn{t) + z{t) (1) 

n—1 

where Soit) and s„(f) are the transmitted symbols at time instant t from the desired user and the interfere!', respectively, with 
average powers i5{|sDp} = Pd and i?{|s„p} = p„. Let VI G C — {1,, L}, denote the Rayleigh distributed random 

envelope of the received signal of the desired user over the diversity path. The output signal-to-interference-plus-noise ratio 
(SINR) r(t) of the system is given by 


r{t) = 


Pollhc 


Pd E ^D.i(0 


N N 

+ E Pn|w(f)^h„(f)|2 TVo + E Pn 

n—1 n—1 


Y 


( 2 ) 


where Y and Z denote the numerator and denominator of the SINR, respectively, w(t) = hD(f)/||hj3(t)|| is the receive filter. 


and Ai „{t) = |w(f)^h„(f)|. Then, A‘^^{t) and Af^{t) are standard exponential random processes with unit means |35 
Given an SINR threshold 7th for the desired user, the level crossing rate (LCR) is defined as 


OO 

LCR(7th) = J 7/r,r (7th, 7)^7, 


(3) 


7=0 


where the time index t is dropped to simplify the notation, and hence, F is the time derivative of the SINR F, and fpp{.,.) 
is the joint probability density function (PDF) of F and F. 

The average duration for which the SINR remains below a threshold, namely the average outage duration (AOD), is defined 
as |34) 


AOD (7,h) 


Fr (7th) 
LCR (7th) ’ 


(4) 


where Fp (7*) = Prob(L < 7th) is the cumulative distribution function (CDF) of F. 

In Sections [111] and [fV] the LCR for the system presented above is derived. For the reader’s convenience, the main derivation 
steps are summarized as follows. 


Main LCR Derivation Steps 

1: The transformations of random variables (RVs) in (|9]l-([TT]i are introduced to find the joint PDF of the SINR and its time 
derivative. 

2; Baye’s rule is used on the conditional PDFs to find the joint PDF of the SINR, its time derivative and the envelopes of 
the interferers’ paths, given in ( [T7] i. 

3: Integration is performed first over the SINR’s time derivative in the LCR integration as shown in ( |T8] l to get ( |2T] i. 

4; The A^-multiple integral in ([2^ is reduced into a double integral through regarding this integral as an expectation of a 
function of exponential RVs as in ( |26| l. 

5; The transformations of RVs are introduced, then their joint PDF is obtained as in ( [30l l to solve ( |26] l. Using 

these transformations, the expectation is rewritten in terms of Qi and Q 2 only as in ( [35] l. 

6; The joint characteristic function of Qi and Q 2 is obtained as in ( |3^ , then their joint PDF is obtained as in ( [34l l. 

7: The joint PDF is substituted in the integral ( [3^ and after a series of binomial expansions and integrations, the final exact 
closed-form LCR expression in (|40li is achieved. 









III. Joint PDF of The SINK and Its Time Derivative 
F ollowing the SINR definition in the time derivative of SINR F is given by 


r = 


2 pl 




( 5 ) 


where „ are i.i.d Gaussian RVs with zero mean. The variances of the latter are tr^ = | |l9[ , where /, „ = ^ 

is the maximum Doppler shift associated with the interfere^ Vn and A are the relative corresponding mobile speed of 
interferer and the carrier wavelength, respectively. Note that for the desired user, is assumed to have the same variance 

<^0 1 = ^ G -C, where ^ is the maximum Doppler shift associated with the desired transmitter (Tx), and 

is its relative mobile speed. This is due to the fact that these variances correspond to the same Tx, and , depends only on 
the relative speed which is the same for all the diversity paths between the same Tx and receiver (Rx). Then, F is Gaussian 


(conditioned on ,, and I, n) with zero mean and variance cr? and its conditional PDF is given as 119 


/ i,Ay '^i,n) 


27r(T2 


V L 


( 6 ) 


where the variance u? is given by 


2 _ 2 I ‘*7 

7 ^2 2 


Ay 


J2. „,2 _2 _ , 472 ^ 


1=1 


2 2 2 
Pnyr.<^n- 


(7) 


It is noticed from 0 that (T? is independent of , and the knowledge of 7 and is sufficient. Hence, 0 is rewritten 


as 


/r|r,A2„(7l7,a,^„) = 


y Ln 


277(7? 


and ( 7 ? is still given by 0 - To obtain fp.t (7,7) in the following variable transformations can be introduced 

F = F 

Un= VneAA 


F = 


Y 


Y 


?Vo + X)n=l No + J2n=lPnUn 


( 8 ) 


(9) 

( 10 ) 

( 11 ) 


The joint PDF can be derived by using the Jacobian approach as follows 

fr,r,Ui,...,UN (7’ 7) Hi, ■ • ■, Un) = ft\r,Ui,...,UN (7l7, Hi,...,«.«) fr\Ui,...,UN ^7,7, Hi, 

where 




(7,7,Hi,...,H„) = -r^ 
I dY 




Ul,...,Ur 


(ui, . . . .Un) = 


1 




I/v ^7 (^No + ^ 


fAy,...,Aj^ (a?, 1 ,-= 


N 


d{Aj„...,Aiy 


/. 


r\r,Ui,...,UN 


(7l7,Hi,...,u„) = 




r|v.(7i.. 


The Jacobians used above are given by 


or 

dr 


= 1 , 


I dr I 
I dY I 


N 

No + X] PnUn 

n—1 


d{Ui,...,UN) 


d{Ai„...,Ajy 


d{Ui,. 

d{Aj„. 

\ 

■ ■,Un) 

-Aiy 

11 

n—1 

n I ; 1 

■ ■ 5 

)■ 

dUi 

dill 

dUi 



■■■ dAfy 

dU2 

dU2 

dU2 



^^i,N 

dUn 

dUt, 

dllN 


dAj, 

■ ■■ 9AI^ 


(13) 


= HI = 1. 


( 14 ) 







































Since the distribution of the sum of i.i.d. standard exponential RVs is known to be an Erlang distribution, then the pdf of 

L 

X = ^ given as 

z=i 


fx{x) = 




rs(i) ’ 

where rs(L) is the Gamma function. Then the pdf of y = p^X defined in (|^ as the numerator of the SINR is given as 


(15) 


T 1 — ^ 

1 

friy) = -rsYjfx {y/Pn) = 


I dY I 

I dX I 


Given the fact that Af^,\/n G N' are i.i.d. standard exponential RVs, and substituting ([^ and into (HI, we get 


(16) 


N 


fr,r,Ui,...,UN (Tj 7) 111) • ■ •) ^n) — 


Xq Priori T ( Xq Pn'^n 


N 


L-1 


T [ No+ Pn 


-T.u„ 


27rCT? 


P^^siL) 


(17) 


IV. LCR Derivation 

By substituting ( [T7| ) in the level crossing rate given a certain SINR threshold 7 th can be obtained as 

Tth-”o ^ 


LCR(7th) = 


7th /■ /■ 

V^Pn^s{L) J J 

UjyJ^—0 

r _-| Tth-^o OO OO 

7,h e f r 


^o + E 


Pn'^n 


n—1 

N 




iVo + E 


Pn'^n 


- E i+- 




y=e dui... dui, 

,=o M 


ay dui... dui, 


(18) 


(19) 


uj^—0 


n—1 


N 


IthNo \ 

L-1 

2 g P£> 

OO 

f 

OO 

f 

N 

i + E 

Pd / 

rs(L 

J • 

■J \ 



un—0 

ui—0 ' 



L —1 iv 




Pn ~ 

Wo^- 



dui 

.. duM 


Pry 1 API(xI 

No NoPd(tI 


^ \ Pd ) 


i—I T-th^o 
e r>n 


v^UL) n (1 + :)^) 

n—1 ^ ' 


( 20 ) 


( 21 ) 


where ( |20| l follows by substituting by (0 in and doing simple algebraic manipulation. By introducing the variable 
transformation «„ = (!+ ^thPu/PiDun, © follows directly where the integral la is defined as 


N 


L-1 


J - N 

un— 0 ui—0 

and the constants a„ and bn are given by 


1 + ^ ^ ClnUn 
\ n=l / 


\ 


N 


1 + E] bnUnS ^"=1 dui . . . dUj, 


( 22 ) 


n—1 


_ Pn/No 

~ 1 I ■yihPn ’ 

Pd 

bn = a„ (1 + £„), and 

^2 


^n — 


^thPn ^ n 

PdoI 


(23) 


Note that from the above definitions, it can be seen that a„,5„ and £„ are all positive Vn G M. 

Now to proceed forward in solving the integral la, new i.i.d. standard exj^onential RVs denoted by Vn, Vn G JV sxe 
introduced. Hence, their joint PDF is given by ... ,Vk) = Define giVi,... ,Vf,) as a function of 

these RVs as 


N 


L-1 


g{Vi,...,V^) = 1 + Eanl4 


\ 


N 


1 + E/ ^rtVi- 


(24) 


























Then the expectation of this function is given by 

oo oo 

E {g (Vi,... ,V„)} = j ... J g{vi,...,Vr,) fvi ... ,t;„) dvi...dv^ 


vn—0 i;i=0 


L-1 


\ 


N 


1 + bnVnC ^"=1 dvi . . . dVi, 


(25) 


^ / N ^ 

= / ... / 1 + ^ CinVn 

vj=0 v,io V n=l / 

where we notice that the integration in ( |25] l is exactly the same as la in ( |22| ). Hence, only for the sake of solving la, we shall 
regard the integration variables Un, Vn G Af as i.i.d. standard exponential RVs and thus la can be rewritten as 

, L-l 


N 


4 = < 1 + 


\ 


N 


1 + bnUn 


In order to find the above expectation, the following RV transformations are introduced 


N 


Ql — ^ ^ b^Un 


n—1 

N 


Q 2 - 51 


an - 


aiK 

bi 


Un, and 


(26) 


(27) 


(28) 


Qn = anUn, n = 3,...,N. (29) 

After finding the domains of the new variables as derived in Appendix 0 the joint PDF /qi,,..,q„( 91 ,... , 9 «) is given as 


/Qi.....Qn('7i> • ■ ■ .9«) = TJ|-/(7i....,[/«(«!, • ■ • ,Mjv) = TT]- 


1 - 57 Un 1 - J2 OnQn 


J 


J 


where | J| is the Jacobian. | J| and are given as 

|J| = 


^(Ql ,...,QAr) 

d{Uu...,UN) 


N 


= (0261 - 0162 ) R an, and 


— 


bi ’ 
bi-b2 


n—3 

n = 1 
n = 2 


a2bi —aib2 ’ 

(b2-bi)(a„bi-aib„) _ b„ , J_ „ _ o at 

anb\{a2b\—aib2) o^nbi ' ’ 5 ■ * • 5 


(30) 


(31) 


(32) 


After tedious algebraic manipulation, detailed in Appendixthe joint characteristic function (Pq^^q^{si, S 2 ) is written 


as 


N N 




1 , 02 ( 51 .S 2 ) = ^ 5^^451 




(33) 


t^n 

where Si, and S 2 are the Laplace variables and fn,gn, 4, At A and ipt n are defined in Appendix [B] by equations ( |9^ , ( |94l i, 
@ and ( |100| i and further simplified in Appendix by equations ( |102| i-( [T07l ). 

In order to find the joint PDF fQi,Q 2 {Qi,g 2 ), double inverse Laplace transform over si and S 2 is performed, resulting in 
the closed-form expression 


N 


/Qi.Q2('?1.92) = a 5^ 4 5Z V't.r 


g-Sn9lg-At,„(g2-A9l) 


(34) 


n—1 t—1 

t^n 


where this expression is valid under condition 92 > /n 9 i,Vn G Af, with emphasis that /i = 0. These conditions result from 
the inverse Laplace transformation and will affect the range upon which Q 2 is integrated as will be seen later. 

The joint PDF of Qi and Q 2 is then used to solve the integration in ( |22| ) as follows. Using the RVs transformation (|27li-([29]l, 
can be rewritten as 


7a = A’ jv^l + 9i (^1 + +72^ I 


(35) 


L-l 

-E 


L-l 


L-l-k 

E 

m—O 


L — 1 — k\ ( ai 
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^-1 
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L — l — k—r 


/+ ( 91 . 72 ) 6(92 rf9i, (36) 
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Ql Q2 




















where ( |3^ follows from a series of binomial expansions of the bracketed term in (|3^. 
Substituting ([3^ in ([3^, Ib can be rewritten as 


where 


h = A 
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N N 


<^1 '>Pt,nl2 

t —2 n —2 t —1 
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under condition that At i > 0 and gn > 0 which are discussed in detail in Appendix D rinc(a, x) = is the upper 

incomplete Gamma function Sec. 8.350, Eq. 2]. We note that the integration limits are derived based on the conditions 
imposed by the double inverse Laplace transform as detailed in Appendix along with the detailed derivation of the above 
integrations. 

Finally, using (|2T]), (|3^, (|J7|l, (|38|l, and ([39ll, the closed-form expression of the LCR can be written as 



where ^ is defined as 
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since rs(L) = {L — 1)1. 


V. Special Cases 


A. Interferers with Equal Powers and Equal Speeds 

In case all the interferers have equal transmitting powers and move with the same speed, i.e. = p, and /n = /i, Vn S JV, 
the LCR can be significantly simplified as follows 
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We note that a„ = a, Z)„ = b, and £„ = e, where they are still defined as in ( |23] l but with replacing p„ and /„ by p, and /j, 
respectively. The detailed derivation of ( |42l i is given in Appendix 





















B. Interferers with Equal Powers and Equal Speeds in an Interference-Limited System 
In an interference-limited system, the LCR expression in (|42li can be further simplified to 


LCR(7th) = 
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where flth — PD/ 7 thPi- The detailed derivation of ( |43l ) from ( |40l i is given in Appendix]^ It can be seen that ( |43] l is in agreement 
with (18) in for the Rayleigh fading scenario. 

Due to the simplicity of the LCR expression of this case, we derive next an expression for the 7 maxj which is dehned as the 
7 th at which the maximum LCR value occurs. This would give us insight on the behavior of the LCR w.r.t. the speeds of the 
desired user and the interferers, i.e. their associated Doppler frequencies. Taking the derivative of ( |43] l w.r.t. 7 *, then equating 
it to zero results in a quadratic equation in 7 * whose solution is 




(44) 


where the negative root was ignored since 7 th has to be positive since it is an SINR value. From ( |44l i we see that as /□ increases, 
7 max decreases, while it increases as /, increases. Moreover, the dependence on the Doppler frequencies is always in the form 
of the ratio /d//i which suggests that if the desired user’s and interferers’ speeds (or equivalently Doppler frequencies) increase 
by the same ratio, 7 niax does not change. This has been confirmed through simulations but the figures are not presented in this 
paper due to space limitation. While ( |44| l is valid for the special case of interferers with equal powers and equal speeds in 
interference-limited systems, nevertheless it gives us insight on its behavior in more general cases as will be shown in sections 
IVII-BI and IVira 


C. Single-Antenna Receiver 

For a receiver with a single antenna, L = 1, (|40|) simplifies to 
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The first equality •ES follows from direct substitution of L = 1 in ( |40l i and the second equality follows from some 
algebraic manipulation as detailed in Appendix [G] Equation ( [46] ) agrees with the LCR expression reported in |19|. 


VI. Approximated LCR Expression 

Eor the sake of simplifying the LCR expression in ( |40l l and getting insight on the behavior of the LCR w.r.t. different 
parameters easily, an approximate LCR is derived in this section. In Q, Eukawa et al. derived LCR expressions for a lower 
bound of the SIR of a point-to-point system subject to multipath with MRC for two different cases. The first case is referred 
to as equal average power (EAP) where the components of the desired signal and interference signals have equal average 
powers, i.e., , = Po, V( S C and p„ = Pi,yn G Af. The second case is referred to as unequal average power (UAP) where 

the components of the desired signal and interference signals have unequal average powers, i.e., p^, 7 ^ j for I f j and 
Pn 7 ^ Pm for n m. The LCR expression is exact for the EAP case while it is an approximation for the UAP case, however 
for both cases LCR expressions are derived for a lower bound of the SIR as mentioned earlier. Next the approach in Q is 
utilized to hnd the approximate LCR of our system. The SINR E in can be rewritten as 
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where R^ and Rf denote the numerator and denominator of the SINR, respectively. Eurther, as shown in ||^, the LCR in ([^ 
can be rewritten as 
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where roj ?"i) ?^i) is the joint PDF of the envelopes of the total desired and total interference signals and their 

time derivatives. It is given by 
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where the first equality follows from the fact that the desired and interfering signals are independent, while the second equality 
follows from the independence of each signal and its time derivative as explained in Appendix along with the reason of 
using the approximation sign. The marginal PDFs in the above equation can be obtained as 
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The PDFs ( [50pT] i can be obtained by properly adapting the variables in equations (23) and (24) of the EAP case obtained 
in 1^ to match the notations in this paper. Note that an important difference between the system in this paper and in Q is 
that here multiple interferers possess multiple maximum Doppler frequencies due to different speeds while in Q there is one 
maximum Doppler frequency only since all components belong to the same interfering user. Another important difference is 
that here a more general case is considered by obtaining the LCR of the SINR while in Q the LCR of the SIR is obtained, 
i.e. the system under consideration in ||^ is an interference-limited system. The details on the derivation of ( [50p^ can be 
found in Appendix [H] 

By substituting ( |49| ) into ( |48l l, it is found that LCR(7th) = AB, where 
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The above derivation follows similar steps as in ||^ Appendix B] and is not included here for the sake of brevity. Subtracting 
from (|56|), A is given as 
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However, for the derivation of B, a different approach from ||^ is needed as shown next. By directly substituting fR,^{ro) 




















and 


in respectively, B can be obtained as 
B = 


2 


N I 

2/i„ei 




Pn 


-TlhPn+Pr) J.2 

r, e popp ' dri 


r\ — y/No 


2v^7, 


L-i AT 


L-i -Wp 


^ePpTinc L + 


1 7ih-/V'o 
2’ PD 


No 

Pn 


^ ' n=l 


(7thPn +Pd) 


L+i 


( 59 ) 


(60) 


where |36 Sec. 3.381, Eq. 9] is used to solve the integral in 

Hence an approximate closed-form expression for the LCR of the SINR given in Q is obtained by the multiplication of 
and ( |60l l. Thus, denoting the approximate LCR given an SINR threshold, 7 th as LCRappiox{ 7 th), it can be written as 
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As mentioned earlier, this approximate expression provides simplicity to get more insight on the behavior of the LCR w.r.t. 
different parameters easily. Lor example it is much easier to plot the behavior of the LCR w.r.t. N and L from the approximate 
expression (|M]l than from the exact one (|40ll. 


VII. Validation of Theoretical Results 

In this section we demonstrate the validity and accuracy of both the exact and approximated LCR expressions derived in this 
paper, ( |40| i and ([6T]i respectively, through comparing them with simulation results for different scenarios. We note that due to 
the complexity of the LCR expression and the fact that it is highly non-linear function in these parameters, and the abundance 
of parameters affecting the LCR and consequently the very wide variety of combinations of these parameters that need to 
be investigated, it takes a lot more examples and complex cross-analysis to investigate the interrelated effects of the different 
parameters of the system on the LCR’s behavior. Since our main focus is to validate the novel theoretical LCR expressions 
presented in this paper and not the cross-analysis of the effect of the system’s parameters on the LCR, we suffice by presenting 
a few examples showing some interesting observations on the effect of the powers and speeds of the desired user and interferers 
on the LCR behavior, while detailed analysis of the effect of different system’s parameters on the LCR is left for future work. 
To show the validity of our LCR expressions, specially the exact one ( |40l l, we simulate three different systems with different 
number of Rx antennas L, interferers N, powers and p^, Doppler frequencies and /„ and SNR p^/No- The values of the 
aforementioned parameters are chosen randomly across a range of possible values in practical cellular system. 

Lor the simulation results, the Rayleigh fading channels are generated using the method described in nzl which is a 
variation of Jake’s method. The simulations are performed using Matlab R2014b, and the results are averaged over 400 channel 
realizations each for time duration of 5 sec with 1 MHz sampling rate which is much larger than the maximum Doppler 
frequencies simulated, hence, capturing the channel’s variations. It is also worth mentioning that since we are performing 
baseband simulations, then the sampling frequency’s value is determined by considering the simulated Doppler frequencies’ 
values only and not the carrier frequency’s. 


A. Effect of Desired User and Interferers Powers 

In Lig. [T] the LCR is plotted for an interference-limited system where L = 2, N = 2, for different values of the interferers’ 
powers, pi, and different values of the desired Tx power po. In particular, five different combinations of interferers’ powers are 
plotted, namely when the interferers’ powers ratios, Tp = Pi, 2 /Pi,i are equal to —10 dB, —3 dB, 0 dB, 3 dB, and 10 dB. Each 
of these cases is simulated for two values of p^, namely 10, and 1, corresponding to dominant and non-dominant desired Tx 
power cases, respectively, w.r.t. the maximum of the interferers powers. In all systems, the maximum Doppler frequencies of 
the interferers, /i „, are 32 and 162 Hz, for n = 1, 2, respectively, while that of the desired user, /□, is set to the minimum of 
both, i.e. 32 Hz. These frequencies correspond to vehicular speeds of 10 and 50 km/hr, respectively, for a cellular system with 
carrier frequency 3.5 GHz, which is band 22 in the LTE LDD bands |38 Table 5.5-1]. 

The LCR is plotted versus the SINR threshold 7 * normalized w.r.t. the average SINR js.vg which is obtained through 
simulations. In the figure, we refer to the LCR obtained by equations ( |40l l and as ‘Exact’ and ‘Approx.’, respectively. 
Lig. 0 shows that the exact LCR matches perfectly the LCR values obtained from simulations. On the other hand, there is a 
deviation between the approximated LCR and the exact LCR. 

We also note that as p^ changes, the average SINR value changes, however, the LCR value does not change for a fixed 
7th/7avg- This can be seen in Lig. where the LCR curves are the same for both dominant and non-dominant desired user 
transmitting power cases, where po = 10 and 1 , respectively. 

Lig. 0 indicates that as the power ratio, Tp, increases, the LCR increases. This can be explained as follows. In general 
as the power of an interferer increases, its effect becomes more pronounce. Since the second interferer has a higher Doppler 
















Fig. 1: The effect of power values on the LCR where L = 2, TV = 2, /□ = 32 Hz, and /, = [32,162] Hz. 


frequency, hence when its power increases, more fluctuation is induced in the resultant received signal leading to more frequent 
crossings around a given threshold. This shows in the cases Tp = —10 dB, —3 dB, and 0 dB, where i is fixed to 1 and p, 2 
increases. In the cases Tp = 0 dB, 3 dB, and 10 dB, still the LCR increases yet it is because the power pi 1 of the slower first 
interfere!' decreases while p, 1 is fixed to 1. This observation can be stated in other words; the LCR increases as the ratio of 
the faster interferer’s power to the slower interferer’s power increases, and vice versa. 

In order to analyze the gap between the approximated and exact LCR and how it is affected by the system parameters, we 
define a metric named maximum relative gap, MRG, as 


MRG = 


max [LCRexact TCRapprox] 
LCRexact, MRG 


(62) 


where LCRexact. MRG is the exact LCR at which the difference between the exact and approximate LCR is maximum. 

Table |I^ lists the values of MRG for the cases plotted in Fig. It can be seen from the table that as the difference in 
the powers between the interferers increases, the MRG increases as well, indicating that the accuracy of the approximation 
declines. This is because the approximated LCR assumes that the time derivative of the interference signal is independent of 
the interference signal itself which is not accurate, unless all of the interferers have the same transmitting powers and speeds. 
Thus, for given speeds of the interferers, as the interferers’ powers get closer in values, the assumption of independence of 
the two signals, and consequently the approximation made within the derivation, become more accurate, and vice versa. This 
shows in the low MRG values that the equal-power case exhibits. 

It is also noticed that for the same interferers’ powers ratio, the MRG values are larger when the faster interferer has larger 
power compared to the opposite case. This is obvious when comparing the MRG value in cases of Tp = —10 dB and 10 
dB together and comparing the cases where Tp = —3 dB and 3 dB together. Again this is attributed to the assumption made 
during the derivation of the approximated LCR, where as the speeds and powers of the interferers get closer in value, i.e. 
more homogeneous interferers, the more accurate the assumption mentioned in the previous paragraph is. One metric of the 
interferers homogeneity is the maximum product of the interferer’s power and speed, maxpi „/i „. Hence it shows that as the 

n 

power-frequency product increases like the case of Tp = 10 dB, the more inhomogeneous the interferers are compared to the 
case of Tp = —10 dB, and hence the greater the MRG and the less accurate the approximated LCR expression is. On another 
note, it is seen that still the approximated LCR provides reasonable accuracy over the relatively low normalized threshold 
range, 7th/7avg < 0 dB, as shown from Fig. noting that the high values of MRG tabulated in Table [H] usually occur at 
relatively high normalized thresholds. 


TABLE II; MRG values for the system in Fig. 


Tp 

-10 dB 

-3 dB 

0 dB 

3 dB 

10 dB 

MRG 

11.4% 

7.3% 

3.7% 

11.5% 

19.9% 


B. Effect of Desired User’s Speed 

In this subsection, the effect of the desired user’s speed of motion, i.e. its maximum Doppler frequency, is investigated. 
Fig. 1^ plots the exact and the simulated LCR for different speeds of the desired user in an L = 2, A^ = 4 system where 
















Po = 10 while the interferers’ powers are p, = [0.07,0.1,0.05,0.12] and the SNR po/No- This setting is chosen such that the 
desired Tx power is dominant w.r.t. the interferers’ to ensure that the effect of the desired user’s speed is evident. On the other 
hand, the choice of the interferers’ powers is irrelevant since they are hxed over the different values of f^, and indeed the 
same behavior was observed when the system was simulated for different combinations of pi. The corresponding noise power, 
No = 0.1, is chosen that way to be comparable to the interferers’ powers to demonstrate the accuracy of the LCR expression 
( |40| i when taking the noise power into consideration. The hgure again shows that the LCR values obtained from equation ( |40| ) 
match the simulations. The curves of the approximated LCR obtained from equation ( |M] l are not plotted for the sake of clarity 
of the figure, however, it is worth mentioning that they are still close to the exact LCR curves. 

The values of plotted correspond to the speed values of approximately 3, 30, 50, and 100 km/hr, hence, ranging from a 
pedestrian Tx to a high-speed moving vehicle. The hgure shows that for a given combination of interferers’ speeds, as the 
desired user’s speed increases, the LCR increases. This is because more Huctuation is introduced to the received signal due to 
the higher Doppler frequency of the desired Tx, causing faster crossings around a given threshold. 



Fig. 2: The effect of fo values on the LCR where L = 2, N = 4, and /, = [30,160,65,100] Hz, po = 10, 
p, = [0.07,0.1,0.05,0.12] and p^/No = 20 dB. 

An interesting hnding is that as the desired user’s Doppler frequency increases, the normalized threshold at which the 
maximum LCR occurs, 7max/7th, decreases, i.e. shifts to the left as Fig. [^indicates. This is shown explicitly from the values 
of 7max/7th mentioned for each curve. This is in accordance with ( |44l ) though this equation is derived for interference-limited 
systems with equal-power and equal-speed interferers. This indicates that although the analytical expression of 7niax is different 
and more complicated in case of systems with signihcant noise power and interferers with unequal powers and speeds, a similar 
behavior is observed as in the simplihed system of (@41, hence it provides a useful insight on the general LCR behavior. It is 
worth mentioning that a high value of LCR indicates that the signal crosses the associated threshold frequently which implies 
that the received signal spends more time in the vicinity. Consequently, the figure shows that as the desired user’s speed 
increases, the received signal spends more time at a lower SINR value. This indicates that the threshold should be lowered as 
fo increases to maintain same outage probability. 

C. Effect of Interferers Speeds 

In this subsection, the LCR is simulated for a system with L = 3 and N — 2 while hxing po = 1, and /□ = 100 Hz. The 
interferers have close powers where p, = [0.1, 0.09], and equal Doppler frequencies /i_i = /,_2 = f, while the SNR Po/No = 10 
dB. The Tx power, po, could have been chosen differently and still the same behavior would be observed as long as it is hxed 
over the different combinations of the interferers’ speeds simulated. Again, the curves of the approximated LCR obtained from 
equation ( |M] l are not plotted for the sake of clarity of the hgure, however, it is worth mentioning that they are still close to 
the exact LCR curves. 

In Fig. the LCR is plotted for fours systems where the speeds of the interferers are given by approximately 3, 30, 50, and 
100 km/hr. Again the hgure shows that as the speeds of the interferers increase, the LCR increases, however, the normalized 
threshold at which the maximum LCR value occurs increases as opposed to the case in the last subsection. Again this is in 
accordance with ( |44| i which was derived for interference-limited equal-power and equal-speed interferers which conhrms that 
this special case is insightful even for cases with medium to low SNR and not interference-limited as in this system ,where 
No = 0 . 1 . 
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Fig. 3: The effect of /, values on the LCR where L = 3, TV = 2, /□ = 100 Hz, /,i = /,2 = f,, Po = 1, p, = [0.1,0.09], and 
p^/No = 10 dB. 


VIII. AOD Derivation 

In this section, the AOD of the system under study is derived. As AOD is defined in Q, the CDF F'r( 7 ) needs to be 
obtained. In p^ , the outage probability Prob{F < 7 th} = Fr( 7 th) was derived for Rayleigh distributed desired signal and 
Nakagami-m distributed interference signals with distinct unequal powers. By substituting m = 1 in [39) , where m is the 
Nakagami parameter, and adapting the notations, the CDF of our system ^^( 7 *) is given as 
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upper incomplete gamma function rin(,(TO + l,a;) [36 Sec. 8.352, Eq. 2] 


where the second equality follows from some algebraic manipulations, and from the following series representation of the 
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By substituting ( |63| l and [40) 1 in Q, an exact closed-form expression for the AOD is achieved. 

Fig. 1^ plots the AOD vs the normalized threshold for the interference-limited systems with {L,N) = (2,2) and {L,N) = 
(4,4). The maximum Doppler frequencies of the interferers are /i = [32,162,65,97] Hz, selected across a range of possible 
Doppler shifts in practical cellular system. Their corresponding transmitting powers are p, = [1,0.5, 0.8,0.3], while for the 
desired user fo = 32 Hz and po = 1. The two-interferer system consists of the first two interferers only. The figure compares the 
AOD obtained from the simulation with the exact and approximate AOD expressions resulting from the exact and approximate 
LCR expressions ( |40| l and •HD, respectively. As expected, the AOD increases monotonically as the threshold increases since 
the probability of the received signal SINR being greater than the threshold decreases as the threshold increases. This results 
in spending more time in fade which in turn increases the AOD. We note that the exact expression matches almost perfectly 
the simulation results whereas the approximate solution is very close to the simulation as well. It is also observed that as L 
and N increase, the AOD increases as well as also observed in ©■ 


IX. Application 1: PER Calculation 

In this section, the PER is derived for uncoded transmission. Ideally, the calculation of PER requires the demodulation of 
all bits in the packet to determine if all bits in the packet are received correctly, else the packet is considered erroneous. 















Fig. 4: The AOD vs normalized threshold for interference-limited systems {L,N) = (2,2), and {L,N) = (4,4), /□ = 32 Hz, 
/i = [32,162,65,97] Hz, = 1, p, = [1,0.5,0.9,0.3], and p^/No = 50 dB. 


However, this mandates a large computational time in a conventional computer simulation. Hence for simplicity, the data 
packet is considered to be received correctly only if the received SINK is above a certain threshold 7 * during the whole data 
packet duration Tp^t Q, as shown in Fig. If at any instant, the SINR drops below the threshold, the received data packet 
is considered to be erroneous. This process can be modeled by using the two-state Markov model Q. The system is in good 
state G if the SINR is greater than or equal to the threshold 7 th, which can identify a target PER. The system is in bad state 
B otherwise. Hence, PER Pg is the probability of being in state B. In other words, Pe = Pb = ^ — Pg, where Pb and Pg are 
the probabilities of being in states B and G, respectively, also known as steady state probabilities. Pg is given as Q 

Pg = ( 66 ) 

where Pcf = 1 — Fr(7th) is the complementary CDF To simplify the notation, define Nc = LCR(7th). Substituting the LCR 
and the CDF derived in previous subsections, the PER of the uncoded system is obtained. It is worth mentioning that the PER 
considers the channel time variations since it is derived using the ESMC model which captures the channel time variations 
and correlation, second-order channel statistics, through LCR. 



Eig. 5: The received SINR over a time interval where L = 2, N = 2, 7avg = 7 * = 15 dB, and = /i_max = 35 Hz. 


In Eig. 1^ the PER plots of two systems with {L,N) = (2,2) and (4,4) are compared at the same normalized threshold 
value of —14 dB as the packet length changes. The PER is plotted versus the packet length normalized w.r.t. the maximum 










Doppler frequency in the system. The desired user in both systems has a transmitting power of pu = 1, and a maximum Doppler 
frequency equal to 32 Hz which is equal to the minimum Doppler frequency among the interferers in both systems. The powers 
of the interferers are in the same order, thus no single interferer is dominant over the others. In particular, the powers of the 
interferers are 1,0.5,0.8, and 0.3 and their corresponding maximum Doppler frequencies are 32 Hz, 162 Hz, 65 Hz and 97 
Hz. The two-interferer system consists of the first two interferers only. We also plot the two scenarios; interference-limited 
system (po/No = 50 dB), and noise-limited system, {p-o/No = —3 dB). For the interference-limited systems. Fig. shows that 
the theoretical PER derived through using the exact LCR expression matches the simulation results very well. It also shows 
that the PER calculated via the approximated LCR expression matches very well, almost perfectly, the actual PER values 
obtained through simulation. This is because the approximate LCR is almost the same as the exact LCR for this system at 
7th/7avg = —14 dB, as shown in Eig.[^ In fact, it is noticed that the approximate LCR is almost identical to the exact LCR over 
relatively low values normalized threshold, which are the typical range of values used in practical systems. The discrepancy 
between the approximate and exact LCR appears most around 7th/7avg = 0 dB, and then diminishes again for relatively high 
values of normalized threshold, i.e. where the LCR value itself decreases again. Thus, for practical systems, the approximate 
LCR expression can be used with high accuracy. The same behavior is experienced for the noise-limited systems. 

Fig. also shows that as the packet length increases the PER increases as well. This is intuitive since the probability 
of having at least one erroneous bit increases. It also shows that as L increases, the PER decreases though the number of 
interferers increases. The PER improvement is attributed to the increased diversity rather than increase in the SINR. Moreover, 
the amount of improvement depends on the speeds, i.e., Doppler frequencies, of the desired user and the interferers as well as 
their relative powers. 



Eig. 6; The PER vs normalized packet length, * /d, for (L, N) = (2,2) and (4,4), where /□ = 32 Hz, /, = [32,162,65,97] 
Hz, /d = max{/D, /i.max}, Pd = 1, Pi = [1, 0.5,0.8, 0.3], and 7th/7avg = -14 dB. 

It is worth noting that this model is very useful to system designers due to the accuracy of PER calculation compared with 
the simulation results. Hence, a large amount of simulation time can be saved, specially for large number of interferers, receive 
antennas, and long packets. We stress that this accurate result is due to the exact LCR expression proposed in this paper. 


X. Application 2; Throughput Maximization for an ARQ-Based System 


A. ARQ Basics 

Many practical wireless packet systems deploy ARQ schemes to improve the transmission reliability. In this section, the 
stop and wait ARQ (SW-ARQ) scheme with unlimited number of retransmissions is deployed in the system under study. In 
such a scheme, during a single ARQ round, the transmitter transmits a packet and waits for an acknowledgement (ACK) from 
the receiver. If the transmitter receives an ACK, then the transmitter transmits the next packet, otherwise it retransmits the last 
packet. The duration of a single ARQ round is given by Tarq = Ts{mi + mo), where Tg is the transmitted symbol duration, 
and mt is the number of symbols per data packet, i.e., data packet length. TOq is the equivalent number of symbols induced 
by the ARQ protocol overhead, e.g., ACK transmission duration and the guard interval for processing. 

Eor an SW-ARQ scheme with unlimited number of retransmissions, assuming a PER Pg and that the receiver detects any 





data packet error that occurs, then the average number of retransmissions M is given as 
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fc=i 


l-P. 


M = (1 - Pe) + 2 Pe(l - Pe) + 3 PI{1 - P,) + ■ ■ ■ + k P^-\l - P,) + ■ ■ ■ = {1 - P,)Y, fc(Pe)""' = 

Then the throughput R of such a system is given by 


rr^ _ TOt (1 - Pe) 


(67) 


( 68 ) 


Plot (wt + mo) Ts ' 

where mr is the total number of user’s information symbols to be transmitted, and Ttot is the total transmission time given by 

mr - 


Ttot = — M Parq- 
rrit 


(69) 


B. Optimal Packet Length 

Since Pg = 1 — Pg, then from (|6^ and (|68|), the throughput can be rewritten as 


R{mi) 


_ I^cTs 

PcF mt e ^cF 

(mt + mo) Ts 


(70) 


where the data packet duration is Tpkt = mtPs- Thus the throughput of the ARQ-based uncoded multi-user system under 
correlated fading channel is a function of the data packet length mt. In the following, the optimal data packet length that 
maximizes the throughput of such a system is obtained as a function of the system parameters by taking the derivative of 
R{mi) w.r.t. mt and equating it to zero 


dR{mi) 

drrii 


g-JV,T,/PcF 

(mt + mof 



N(.moral 1 = 0 . 


By solving the quadratic equation (|7T]l, the optimal data packet length mopt is given as 


(71) 


mo 

mogl — 



4Pcf 

moN(,Ts 



(72) 


where Pcf and A), are functions of 7 th and the system parameters, i.e., the powers of the desired and interfering users, and 
Pn, respectively, the maximum Doppler frequencies of the desired and interfering users, /□ and respectively, the number 

of interferers N, and the number of diversity paths L. 

We note that while in this paper we are considering obtaining the optimum packet length for one receiver, the same framework 
can be used in broader networks such as interference networks, where each transmitter has only one intended receiver. This 
can be done by obtaining the optimum packet length for every receiver in the network through following the same approach 
presented here for one receiver. Note that the packet length will not be the same for each receiver, since for each receiver, there 
is a different desired transmitter and the transmitters in the network have unequal powers and speeds. Thus, for every receiver, 
the corresponding powers and speeds of the desired transmitter and interferers are different, hence the optimum packet length 
differs for every transmitter and receiver pair. 


C. Discussion 

From ( |70| , it is seen that the data packet length mt has adverse effects on the throughput. To elaborate, the throughput is 
rewritten as R(mt) = =^/i(mt)/ 2 (mt), where 

/i(mt) = —^—, and / 2 (mt) = (73) 

mt + mo 

It can be seen that /i(mt) is monotonically increasing while / 2 (mt) is monotonically decreasing with mt. However, the 
steepness of each, i.e. how fast the curve increases or decreases, depends on the values of mo and N(.Ts/Po^, respectively. 
As mo increases, the slope (steepness) of /i(mt) decreases, thus taking longer time to reach the asymptotic value 1. This is 
confirmed by (0, where for a fixed N(.Ts/Pcf, as mo increases, the term inside the parenthesis decreases while the term 
outside the parenthesis increases. However, the effect of the latter is greater, thus resulting in increasing mopt as mo increases. 
This result agrees with the intuition to increase the data packet length as the overhead increases. Note that this result does 
not contradict with the fact that as the data packet length increases, the packet bears more channel changes thus increasing 
the probability of error. This demonstrates the importance of obtaining an optimum data packet length that is long enough to 
overcome the large overhead and small enough to reduce the PER. On the other hand, as mo decreases, mopt decreases as 
well, until the point where the throughput is monotonically decreasing, depending on the relative values of mo and NcTs/Pcf- 
In this case, the optimum data packet length would be the minimum data packet length allowed by the system. 












Fig. 7; The normalized throughput vs packet length in symbols for {L,N) = (2,2) and (4,4) interference-limited systems 
where p^/No = 50 dB, /o = 32 Hz, /, = [32,162,65,97] Hz, = 1, P: = [1,0.5, 0.8,0.3], and 7th/7avg = -14 dB. 


With all other parameters fixed, as 7 th increases, the complementary CDF Pcf increases monotonically while Nc increases 
until around the average SINK value of the received signal then starts decreasing as 7 th further increases as shown in Fig. 
The relation between Wopt and 7 th is not straightforward, which needs further investigation through ( |72l i. 

Fig. I?] plots the throughput normalized to the data rate (Rs = 1/Ts), i.e. R{m)/Rs, for two interference-limited systems, 
namely the two- and four-interferer systems described in Section for two data rates of values 10® samples/sec (100 ksps) 
and 10® samples/sec (1 Msps) at normalized threshold value of —14 dB. In these systems, the ARQ protocol overhead is set 
to TTio = 100 symbols. The figure shows that the exact and approximate theoretical throughputs, obtained via the exact and 
approximate LCR expressions and the FSMC model, match very well the simulation results. 

We note that although the throughput is normalized in the figures, as the data rate increases, the normalized throughput 
increases. As seen in ( fTOl i, the throughput is dependent on the data rate, Rs = l/P, in the exponential function that is a resultant 
term from the PER expression. Thus, as Rg increases, Tg decreases and the normalized throughput increases as obvious from 
( |70l i given a fixed number of symbols in the packet and a fixed SINR threshold. 

As for the four-interferer system, it is noticed that the normalized throughput is very close for both data rates and the 
approximation is pretty close to the exact values. This is because as noticed from Fig. th e PER for the (4,4) system is 
small compared to that of the (2, 2) system. Hence the 1 — Pe factor in the throughput in^681l approaches 1 and the effect of 
the data rate, which is inherent in the Pg expression, is diminished in the normalized throughput. 

Fig. 0 also shows that there is indeed an optimum packet length TOopt which is more evident in the two-interferer system 
curves. In particular, rriopt obtained from (0 for the two-interferer system with data rates 100 ksps and 1 Msps by using the 
exact LCR expression are 1, 328 and 4,304 symbols, respectively. It is obvious from the figure that they match the simulation 
results very well. It is intuitive that as Rg increases while fixing all other parameters, mopt increases as well which agrees with 
This is because fixed Ng and Pcf correspond to fixed Doppler frequencies in which case as the symbol time decreases, 
more symbols can be included within one data packet such that the data packet time duration Pp^t and the normalized data 
packet length Ppkt/d are still same, where fd = max{fo, Thus the packet experiences the same fading behavior as 

with larger Tg and less number of symbols. It is worth mentioning that this argument implies linear increase of mopt w.r.t. Rg, 
however, this is true only if the ARQ overhead mo is increased by the same factor by which Rg increased. Recall that mo 
incorporates the ACK transmission duration which is dependent on the number of symbols in an ACK packet which can be 
fixed and is independent of the data rate. Thus, mo cannot increase linearly with the increase of Rg and hence, neither can 
mopt. This is clear when examining the effects of mo and Tg on mopt in ( |7^ . It also shows in the mopt values obtained for the 
two-interferer system in Eig. where for a fixed mo = 100, when Rg increases by a factor of 10 from 100 ksps to 1 Msps, 
mopt increases from 1, 328 to only 4, 304 and not to 13, 280. 

Eor the (4,4) system, the normalized throughput seems to keep increasing and not decreasing afterwards, however, this is only 
because of the short packet lengths over which the simulation was performed. Indeed, the theoretical optimum packet lengths 
for this system with data rates 100 ksps and 1 Msps are obtained from ( [72l i to be 6,450 and 20, 503 symbols, respectively, 
which are beyond the simulated range in the figure. 








XL Conclusions 


In this paper, the exact and approximate closed-form expressions of the LCR and AOD are derived for multi-user 
multi-antenna wireless communication systems under AWGN and fading channels. The effects of system parameters, i.e., 
the desired user’s and interferers’ powers, and their speeds, on the LCR of the system are analyzed. The derived expressions 
are used in two different applications via FSMC model, evaluating the PER, and maximizing the throughput for an SW-ARW 
based system by setting the optimum packet length as a function of the system parameters. The theoretical expressions are 
validated through comparisons with the simulation results. Extending the model presented in this work to coded transmission 
systems, and throughput maximization via joint optimization of the packet length and modulation size are deferred for future 
work. 


Appendix A 

Solving the integration 


In this appendix, the detailed derivation steps to solve the integral in ( [2^ are provided. As mentioned in Section IV 
to solve the integral in ( [22l i, t/„,Vn G Af aie assumed to be i.i.d. standard exponential RVs, then their joint PDE is given 
by ... ,Un) = Let p({7i,..., t/„) be a function of these RVs, then its expected value is defined as 

E{g(Ui ,[/„)} =/.../ g{ui,... , ^^dui... du^. Using this argument, la can be regarded as in 
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In order to find the above expectation, the following RV transformations are introduced. 
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The inverse transformation, i.e. in terms of Qn, can be obtained as follows. Prom 


Qnictn, n = 3,..., in (28 1 , U 2 can be written as 
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by substituting for {7„ = 
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Using ( |27] i, ( |29] l and ( fMl i, Ui is given 
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Note that although Un,yn G Af aie defined over the non-negative domain only, i.e. Un > 0, it might seem straightforward 
that Qn > 0, Vn, as they are either just weighted sum of (7„, n = 1,2 or a scaled version of Un,n = 3,N. However, this 
is not always the case for Qi and Q 2 as it turns out. To find the domain of Qi, we substitute for Ui by (f75|) in the inequality 


N 

Ui > 0 , then it follows that Qi > PnQn, where /3„ is defined as 
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and /3i is not defined. 
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Similarly, from ( fMl i it follows that Q 2 > bnQn, where hn are defined 
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while hi and ft ,2 are undefined. One very important note here is that the above domain of <52 is correct only if 0261 > 0162 else 
the inequality is reversed. Por the rest, of the paper this condition is assumed to be always satisfied and in Appendix jP] we show 
the implication of this and how to satisfy it physically given a real cellular system. Pinally, the domain of Qn, n = 3,..., N 
is simply Qn > 0 since they have one-to-one relation with u„. 
















Thus, the joint PDF (gi, • • •, ?«) is given as 
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where | J| is the Jacobian. |J| and can be obtained to be 
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Using the characteristic function approach, detailed in Appendix [b] the joint PDF /qi,Q 2 ( 9 ij 92 ) is given by 
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where 52 > /n9i for every n S Af, and p„, /„, At.„, ijjt,n and A are defined in ( |103| l-( |T07] i. 

Using the RVs transformation ([27]i-(|29]l, (|26ll can be rewritten as 
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where ( |79| l follows from the binomial expansion of the two terms in square brackets in ( |78l l. Similar argument leads to ( |80l l. 
Substituting ([3^ in ([3^, /{, can be rewritten as 
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The integration limits are derived based on the conditions imposed by the double inverse Laplace transform as mentioned 
previously. Note that in ( [JT] ) the term corresponding to n = 1 is separated from the rest of the terms because it is a special 
case where /i = 0 which affects the domain on which <72 is integrated as shown in I 12 and I 21 in ( [ST] ), leading to two different 
integration results as will be shown next. 

It is obvious that Ii can be easily solved since it is decoupled into the product of two integrals In and Hence, from 
136 Sec. 3.382, Eq. 4] and Sec. 3.381, Eq. 4], it can be shown that 


Ill = 


e^^Einc (fc + |,pi) 


fc+l 

9i 


and 


/12 = 


m! 

\m+l ’ 


(82) 


under condition that > 0 and At_i > 0 which are discussed in detail in Appendix [D| 

On the other hand, I 2 is nested, so we integrate over q 2 first, then over <71 since the integration limits of (72 are function of 
<71 as imposed by the inverse Laplace transform conditions. Erom |36 Sec. 3.381, Eq. 3], I 21 and I 2 can be obtained as 
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where the second equality results from expanding the incomplete gamma function in the series form using |36 Sec. 8.352, 
Eq. 2]. Substituting by (|83]) in ([8T]i, I 2 can be solved as follows 
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where the variable transformation y = 1 + gi is used in ( [85] l, while ( | 86 l l follows from applying the binomial expansion on 
{y — 1)’'. Equation ( |39l l follows from solving the integration in ( | 86 l l using Sec. 3.381, Eq. 3], under condition that 9n ^ 0 - 
Einally, using (|2T]), p 6 |, p7li,p8]l and p9ll, the closed-form expression of the LCR can be written as 
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since rs(L) = {L — 1)!. 


Appendix B 

Derivation of the Joint PDE fQi , Q 2 { qi , q 2 ) 

Eirst, the joint characteristic function of Qi and Q 2 , (si, 52 )^ is obtained, where si and S 2 are the Laplace variables. 

To be able to use the Laplace Transform concept, we change slightly the definition of the characteristic function to be 
(/?Qi (si, S 2 ) = instead of the conventional definition (si, S 2 ) = Of course the 





























effect of introducing the negative sign has been taken into account through the derivation by adjusting all the conditions and 
domains accordingly. Using the joint PDF (9i) ■ ■ • >'Zn) in <|30l), can be obtained as follows 


(87) 


<Fqi,Q 2 ('Si,S 2 ) = y ...y e , q„)dqi ... dq„ 


Qn Qi 


1 

\J\ 


pij 

J ... J 


Qn—O 


91= E /5n9r. 


( 88 ) 


1 


(si + ai) I J| 


qn—O 


92= E h„gr, 

71=2 


1 


I J| (si + ai) ((si + ai) /?2 + S 2 + 012 ) 


g-[(lll+ai)/32 + S2+«2]92 g EJ(«l+“l)/3.+«r.]9r. 


00 00 N 

r f - E ((^n/32+/3n)(si+ai)+/t„(s2+Q!2)+a7i)9n 


5jv=0 93=0 

1 


(89) 

dq^... dqtf (90) 
-. (91) 


N 

I J| (si + ai) ((si + Qfi) /32 + S 2 + OL2) {{hnfd2 + Pn) Si + {hn (32 + /3n) CTl + hnS2 + ZlnCr 2 + Oln) 

n—Z 

After obtaining the joint characteristic function of Qi and Q 2 , partial fraction is done over si. Thus it can be rewritten as 
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Note that At^„ are undefined for t = n. Noting that Cn(s 2 ) is a function of S 2 , then it is expressed in partial fraction form 
over S 2 as follows: 


N 


Cn(s2) = 






n t -1- \ t tbi '^2 + At,„ ’ 

n (s2 + At,„) ‘^1 

t^n 


(98) 


where the partial fraction coefficients „ are given by 
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By substituting (|98|) into (|92ll, the characteristic function can be rewritten as 
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In order to find the joint PDF /qi,Q 2 (' 71 j(Z 2 ), the double inverse Laplace transform is performed over si and S 2 - The inverse 
Laplace transform is taken over si first, followed by S 2 as follows: 
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where this is valid under condition <72 > friQi for Vn G Af. This condition results from the inverse Laplace transformation and 
it affects the range upon which 52 is integrated as was seen in Appendix 


Appendix C 
Simplifications 

After some tedious algebraic manipulations, the parameters of (|40ll can be simplified in terms of a„ and e„ as follows: 
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Appendix D 

Notes on the LCR Derivation 

It is worth noting that the LCR expression in ( |40l i is only valid if the conditions mentioned within the derivation are satisfied. 
For the sake of completeness, these conditions are discussed here in detail. The conditions are given as follows: 
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where the condition ( |109| l is essential to have the domain of g 2 > Pn<ln, else the inequality is reversed as mentioned 

_ __ n—3 

earlier. In order to find this domain U 2 is substituted for by ( |7^ in the inequality U 2 > 0 and then both sides are multiplied 
by (a 2 &i — 0162 ) which implicitly means that this quantity is assumed positive. Using ( |2^ , we can write aibi — 0162 = 
0102 (ei — £ 2 ) > 0. Since o„ > 0,Vn, then this condition simplifies to 


£i> £2 ^ Picrf > F 2 CT 2 : 


( 111 ) 


where ( [23] l is used to reach the equivalence on the R.H.S. We note that this condition does not impose any constraints on the 
practical systems that this LCR expression can be on. In fact, this condition can be easily satisfied, without loss of generality, 
by ordering the interferers such that the first interferer has a power-variance product (p„tT^) greater than that of the second 
interferer. 

Conditions ( |108| l arises from ( |8T] i and ( [86] l, since these solutions are only possible under these conditions Sec. 3.382, 
Eq. 4], Sec. 3.381, Eq. 3]. It can be shown that these conditions will always be satisfied. As shown in Appendix 
is simplified as = l/ 6 „,Vn G M and since > 0,Vn G M, hence, is always positive and condition ( |108[ ) is always 
satisfied. 

Einally, condition ( |110[ l is necessary to solve p^ Sec. 3.381, Eq. 4]. Unfortunately this condition is not guaranteed 
to be satisfied for all values of the systems’ parameters. Hence, in the following this condition is investigated along with the 
range of 7* over which this condition is satisfied and thus over which the LCR expression ( |40l l is valid. Erom ( |105| l. A* 1 can 
be written as 

bi — bt ai(l + £1) — at(l + £t) 


At,i — 


t = 2,...,N. 


( 112 ) 


atbi - aibt aiat{£i - £t) 

Thus it is obvious that to satisfy ( |110| i, both the numerator and denominator of ( |112| i have to be simultaneously positive or 
simultaneously negative. Since for t = 2, the denominator has to be positive to satisfy ( |108| ) as well, the interferers shall be 
ordered such that the first interferer has a power-variance product (p„cr^) greater than that of all of the other interferers or 


PitTi > pta^ 


El > Et ^ atbi — aibt >0, t = 2,..., N. 


(113) 


Accordingly, ( |1 10| l is satisfied if the numerator is positive as well, i.e. if 61 > bt,t = 2,..., N. This inequality can be expanded 
in terms of the systems’ parameters using (| 2 ^ and some algebraic manipulation as follows 


bi > bt 


Ppgp + 


Pd 

Pi 


+ 7th 


> 


FpO-d + IthPtcrf 


Pd 

Pt 


+ 7th 


t = 2,...,N. 


(114) 


Erom ( |113[ l and ( |114| i, it follows that ( |110| i is satisfied for all values of 7 th iff Pi > Pt: ^ = 2,..., N. However, there might 
be some interferers in the system where pi < pt for some value of t ^ 1. In this case, we derive the range of 7 th over which 
condition (|110|i, and consequently the LCR expression (|40li, are valid as follows. 
























Condition dl 141) can be rewritten as 


ithPiPt {picrl - Ptf^t) + TthPo {pWi - Pl<^t) + PdC^o {Pi - Pt) > 0. 
Solving the quadratic equation on the left-hand side of (|115[), the roots of 7 th are given as 


7th = 


-Pd {pIoI - Pt(Tt) ± \/Pd (Plf^l - Pt^tf - 4piPtPDO-2 (picrj - ptaf) {pi - Pt) 


(115) 


(116) 


‘^PiPt {pi(tI - ptcr'^) 

Given that the interferers are ordered such that ( |113[ ) is satisfied, it can be seen from ( |116| l that if pi > pt, then the roots 
are always negative and since the roots represent SINR threshold which has to be a positive value (on linear scale), then this 
means that the inequality ( |114| i is satisfied for any positive value of 7 th. However, if pi < pt, then the root with the positive 
square root is positive which means that to satisfy the inequality ( 1 14| i, 7 th has to be greater than this root in value. 

In summary, given that the interferers are ordered to satisfy ( |113| l, the LCR expression ( |40| ) is valid for all values of 7 th if 
Pi > pt,t = 2,..., N. If for Vt G T, Pi < Pt, the LCR expression (|40li is valid over the range of SINR threshold given by 


7th > max 

tGT 


-Pd (p?cr? - p^^a^) + pI (pIuI - pjalf - Apiptp^ul (picr? - ptof) (pi - pt) 


2piPt (Picr? -PtCTt) 


(117) 


On another note, taking a look at the procedure of the derivation and the conditions under which the LCR expression is 
derived, it shows that the LCR expression in ( |40l i is in general valid for equal, unequal power co-channel interferes with 
equal/unequal speeds, and any combination of these scenarios, except the case: interferers of equal powers with equal speeds. 


where this case is derived as a special case in Section V-A 


It is worth noting that the derived LCR expression does not ignore the noise in the SINR definition unlike many other works 
in literature that are valid for high SNR regime or interference-limited scenarios only. In conclusion, the LCR expression given 
in ( |40l i is valid for any combination of interferers powers and speeds, in interference-limited, noise-limited systems and all 
SNR regimes. 


Appendix E 

LCR Derivation eor Intereerers with Equal Powers and Equal Speeds 

In this section, the LCR expression in ( |42l i is derived. Eollowing the same steps as in the general case, the LCR is given by 
(IZB and the integral in ( |26l t can be rewritten as 

L-l 


N 


Ia = E\ l + aJI^ 


\ 


N 


i + bj2'‘ 


(118) 


where a„ = a, = b, and £n = £,^n G JV and they are defined as in (| 2 ^ and still [/„ are considered i.i.d. standard 


N 

exponential RVs. The variable transformation Q = ^ C/„ is introduced, which is an Erlang distributed variable whose PDE 

n—1 

is given as 


/q(7) = 


qN-^e-q 
rs(iv) ■ 


Then the integration in ( |1 18 I 1 can be written as 

I a = E ^{1 + aq) Vl + 

00 

= y (l-f aq)^~^ 7/1 -f bq 


q^-^e-q 


dq 


q=0 


1 


UN) 

J {£ + [I + bq])^~^ Vl + bq q^-^e~q dq 


T,{N){1 + £ 

1 

r,(iV)(i + £)^-i ^ V I 


q=0 

L-l 


^ E J (1 + bq)U e-q dq 


q=0 


eb 


^ ^ / T 1 \ 

E l ^ ~ ^ \ 


r,(7V)(i + £)i-i ^ V I 


(p-ir” y' 




dy 


(119) 

(120) 
(121) 

( 122 ) 

(123) 

(124) 
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eb 


L-l 


r,(N)(l + e)^- 


eb 


L-l N-1 

3iEE 


L-l 


r s 


N-1 


N - 1 


rs(7V)(i + £)^-i ^V i 

1 l^Q m=0 ^ 


m—0 

L-l\ /N -1 


(- 1 ) 


N-m-l L-l-l / m+i+i ly 


^m+<+5g-5y 


!/=l 


TO 


^L-l-1 ("^ + ^ + 2’ b) 




— m—l— 4 


(125) 

(126) 


where \122) follows from ( |121| l by substituting for a = &/(l + e), and ( |123| l follows from the binomial expansion of 
(e + [1 + bq])^~^ in \ \22) and interchanging the order of the summation and integration. The variable transformation y = 1+bq 
is used to get \\2A) from which ( |125| ) follows by another binomial expansion of (y — 1)^”^. Finally, the integration in ( |125| ) 
is given by Sec. 3.381, Eq. 3] resulting in ( |126| l. 

By substituting \\26\ in (|2T]i, the LCR for the equal-power equal-speed interferers case is given by 


LCR(7th) = 


EE 


(- 1 ) 


N-m-l ^L-l-1 


Einc (w + I + f, I) 


\/^(l + s)^ ^ 


n to! {L-i-iy. (N-m-iy. 


— m—L— 7 


( |4^ revisited) 


where the identity rs(a;) = {x— 1)1 for an integer x is used. It is worth mentioning that this expression has never been reported 
in literature up to our knowledge. 


Appendix F 

LCR Derivation eor Intereerers with Equal Powers and Equal Speeds in an Interference-Limited 

System 

To find the LCR in an interference-limited system where all interferers have the same transmitting power and move with the 
same speed, the limit of ( |42| l as A'o —0 can be easily taken. First, ( |42| l is rewritten in terms of No- From ( |2^ , b is rewritten 
as 


1 


b= — 


Pi 


1 + — 
° PD 


1-f 


TthPigf 




(127) 


where the dependence of b on No is explicitly shown. We assert that neither a, nor e is a function of No per their definitions 
in \\21) and (|2^, respectively. Then, the LCR in (|42ll can be rewritten as an explicit function of No as follows 


’ 2 _ 2lh 


LCR(7th) = 


g PD g a(l + e) L—lA^—1 

EE 


(- 1 ) 


N-m-l ^L-l-1 


Line -f ^ -f 2 1 a(i+e)) i^o) 


N-\-L—m—l — 2 


(d (1 -I- e)) 


N-m-l-% 


(128) 


lies within the term A 


, ^r 1 Z! to! (L-Z-1)! (W-TO-l)! 

V^(l -I-( 1 -I- ^ 1 1=0 m=0 ^ ^ ^ ’ 

From ( |128| l, it is noticed that the dependence on A 
e( PD ''’“U+oE°(A^q)^ Tine ^TO + I + ^, where z = N + L — m — I — 2, I = 0,... ,L — 1 and to = 0,..., A — 1. 

It can be seen that z = 0, when (I, to) = {L—, N — 1) and z > 0 for all other combinations of I and to, which leads to 

limA = |“7 3, V.Nl^iL-l.N-l) 

-f Z-f I) , (Z,to) = (L - 1, A - 1), 

where the identity rinc(a;,0) = rs(a;) is used. Then, the LCR of the interference-limited system under consideration can be 
obtained by taking the limit of (|128|l as No tends to zero as follows 


LCR = lim LCR 7,h = -- ^ - ..j/ . . - . 

y^{l I c)^~^ ^1 I (^ ~ !)• (A — 1)! a-i-Ts (1-p e) ■'■2 




1 T^{N + L - \) (7thPi/pD)^ ^ 

r,(A)r,(L) + 

r.(A + £-2) / fyy 
r,(A)r,(L) + 


(130) 


(131) 


( [43] revisited) 


where flth = Po/TthPi- ( |131| l follows from direct substitution for the variables d, and e by their definitions in ( |127| ) and ( [2^, 
respectively. To get to the final expression in (|43]), we also used = tt^/^ and af = which are defined in Section III 



















































Appendix G 

LCR Derivation for Single-Antenna Receiver 

By directly substituting for L = 1 in ( |40l l, the summations over k, m, r and w disappears and we substitute for k = m = 
r = w = 0 and „i = 1, then (|40ll is given as 


Tlh^o 


/2o'o7thAo e ™ A ^^<5 ^ { 2 '’9n) 


LCR(7th) = 


ttPd 


N , 

n 1 

n^l ^ 


1 + ^ 


t^n 


9n ^t,r, 


( 1^ revisited) 


By re-ordering the terms, the LCR can be expressed as 


Tlh-^O 


LCR(7th) = 


12a^'ya^No e pd' 


N 


ttPd 


AT . 

n 


1 1 7lhP77 

PD 


9n^t,n \/9n 

t^n 

^ ^ ^ 

^MERL 


inc J 2 ’ 9n I 5 


( 1^ revisited) 


where this equation is the same as equation (26) reported in | [T^ where denotes Sn in [p~^ to avoid confusion, and g„ 

is equivalent to >V„ in ||l9) which is proven as follows. First, a„ and £„ are substituted for by (|23|l in ( |103| l resulting in 


1-f 


7lhPr» 


gn = 


a„(l-fe„) PiL _L 


= Wn 


(132) 


where the last equality follows from the definition of >V„ in equation (24) in | [T9t . Thus, it is obvious that gn = Wn- 

Second we prove that is equal to the summation in ( |46l l as follows. Using equation (27) in p^ , and the fact that 

gn = Wn, is rewritten as 
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n 


Wt 


N 

n 


9t 


= 1 t=l 9t - gn 

t^n 


N N / \ i\ i\ / \ 

n fx n (i - i) n i>. n (i - i) 

^=1 ^ / t—1 t—1 ^ ' 


N N 


t 

t^n 
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n b, 


t— 

t^n t^n 
N 


t^n 


N 


, „ ^ , b^-^Ubt n (atbn-anbt) 

A A t —1 t — l , t^n 

~ ~Q N ~ ~Q N N 

n {bn- h) ” n {atbn - anbt) n {bn “ &t) 

t—l,t^n t—l,t^n 

=n 

t^n 

From the definition of the partial fraction coefficient Cn(s 2 ) in (|98l), by substituting for S 2 = 0, we find 


(133) 


N 


Cn(0) = 


Combining (|133|l and (|134|i, it can be proved that 
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n 

t^n 




t^n 


'^t,n 

^t.n 


(134) 


t—1 9n^t,n 


(135) 


t^n 


and hence, ( |46] l is the same as equation (26) in p9) which validates our general LCR expression in ( |40l i. 


Appendix H 

Derivations of (|^ - (|^ 

The joint PDF / 4 |^,(r„P,) can be written as = fR,\R,{'bi\ri)fRi{r,), where /fl,(r,) and are derived 

as follows. 

Derivation of fR^irf)'. Let H denote the interference component of the signal, then H is defined as 
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H =^Pn 


n—1 


(136) 




























( 137 ) 


The PDF of H is given in ||^ Eq. 68] for the UAP case, which is the case under consideration in this paper. 


as 


N 


Pn 


77 > 0 


where are constants given by (|54ll. Since i?i = ^/Ng + H, then the PDF fR^r,) can be given 


as 


= fniv) 


dH 

dR] 


N 
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Pn 


n > \/^o- 


N 


Derivation of fR^irf)'- Since i?, = uNg + ^ Pn '^hen its time derivative yields 


([5^ revisited) 


R^ = 


N 

n—1 

V n=l 


V Pn^l,n^l,r, 
n—1 

R> 


Since, as mentioned in Section 


m 


(138) 


Ai^n is a zero mean Gaussian random process with variance af^ = 119 , then 

it follows from ( |138| l that f, is a linear combination of Gaussian RVs ai^n given a, Hence r, is conditionally (Russian 
distributed with zero mean and variance given by 


E{Rn = 


N 

E 

n—1 


E 


N 


(139) 


No + E P™ 


n—1 


It can be seen from the above equation that the variance E{Ftf} depends on and consequently i?, which makes finding 
/^il^i(f,|ri) a complicated process. Instead, the variance is approximated as follows 
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E{Rf} 


E plE{AUal 

n—1 
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E p>l 

n—1 


E /E pI 

n—1 


Ng + E P’T’ E{Af„} Ng + E P" Ng + E 


(140) 


where the fact that Af^ is a standard exponential RV with unit variance is used. Using this approximation, E{Bf} becomes 
independent of i?,, it can be seen that R, is independent of i?i, and hence, the conditional PDF can be written as 




y/2Tr&f 


([53] revisited) 


Remark For it can be easily shown that Rg is independent of Rg via following the same steps but without 

the need of doing any approximations. This is because the powers of all of the components of the desired signal are equal. 
Thus, the PDFs fRg{rg) and fp^{rg) can be easily obtained to be as in the EAP case in equations (23) and (24) in |j^ 

2 -lA 

fRgiro) = — - — ° r e , and 


fAgi^o) — 


(L-1)! pL 
1 






( |50| revisited) 
([5T| revisited) 


III 


where cTq = 2'k‘^ ff is the variance of A^ , for any I G C as previously mentioned in Section 

The accuracy of approximation in ( |53| l is evaluated in Fig. [8]|^ Fig. |8]|^ plot the approximated PDF and the distribution 
obtained from simulation for a system with L = 2 and N — 2 where the two interferers posses maximum Doppler frequencies 
of values 108 Hz and 540 Hz. Defining the ratio between the transmitting powers of the two interferers as ^ = 10 log(p 2 /pi), 
and fixing the power of the second interfere!' P 2 = 1 for the sake of comparison, the simulation was done for the cases when 
f = 0 dB, i.e. equal powers, and ^ = 10 dB for two different values of noise power Ng, namely Ng = 10“^, and 1, representing 
interference-limited and noise-limited systems. The transmitting power of the desired user is fixed to = 1 and consider that 
it moves with the maximum speed of the interferers, i.e. fg = 540 Hz. When the theoretical approximated PDFs ( |5^ for both 
cases ^ = 0 dB and ^ = 10 were calculated, it was found that they have very close values, almost the same, so we suffice by 



























Fig. 8; Comparing approximated and actual PDF of in the interference-limited L = 2, N = 2 system for the cases of equal 
and unequal powers where No = 10“®, pi a = 1, Po = 1, /i = [108, 540] Hz, and = 540 Hz. 


plotting only one of them and denote it by Theory in the legend to avoid crowding the hgure with many curves with different 
markers. When both PDFs were plotted, they were almost perfectly on top of each other, thus it does not matter which one 
is chosen to plot. Fig. shows that the difference between the approximated PDF and the actual distribution is negligible in 
case of ^ = 10 dB. Surprisingly, when = 0 dB the approximated PDF is less accurate. It is worth mentioning here that even 
when the interferers powers are equal the PDF ( |5^ is still an approximation to the exact fp^{ri) unlike what is reported in H 
because the interferers have different Doppler frequencies. Moreover, even if the interferers move with equal speeds, the PDF 
is still an approximation because we consider the SINR while Q considers the SIR only, hence ignoring the effect of noise. 



Fig. 9: Comparing approximated and actual PDF of rj in the noise-limited L = 2, N = 2 system for the cases of equal and 
unequal powers where Ng = 1, Pi a = 1, Pd = 1. /i = [108,540] Hz, and = 540 Hz. 

As the noise power increases, the approximate PDF deviates noticeably from the exact PDF whether the interferers’ powers 
are equal or not. This is obvious in Fig. where No = 1 which is equal to the dominant interferer’s power. Notice that although 
the deviation is noticeable in the figure, the approximation is still good as was shown in the numerical results. 
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